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We present a quantitative investigation of the impact of tip radius as well as sample type and 
thickness on the lateral resolution in piezoresponse force microscopy (PFM) investigating bulk single 
crystals. The observed linear dependence of the width of the domain wall on the tip radius as well 
as the independence of the lateral resolution on the specific crystal-type are validated by a simple 
theoretical model. Using a Ti-Pt-coated tip with a nominal radius of 15 nm the so far highest lateral 
resolution in bulk crystals of only 17 nm was obtained. 
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Ferroelectrics attract increasing attention due to their 
applicability, e.g., for electrically controlled optical ele- 
ments for efficient frequency-doubling 0, [If, for pho- 
tonic crystals [H, or for non- volatile memories with an 
otherwise un-reached data-storage density The size 
of the domain structures required for those applications 
varies from a few microns down to some nanometers. 
The smaller the domains are the more the properties 
of the domain walls become important. Generally the 
width of the domain walls is expected to be a few crys- 
tal lattice cells @, 0, H|- Measuring the birefringence 
with scanning near-field microscopy in LiTaC>3, however, 
showed a distortion of the crystal over a width of 3 pm 
across the domain boundary 0. Detecting the Raman 
modes with a confocal defect-luminescence microscope, 
an influence of the domain wall in LiNbC>3 further than 
1 pm into the surrounding material was observed [l(| • 
On the other hand, high-resolution transmission electron 
microscopy yielded a domain wall width in PbTiC>3 of 
< 2.5 nm and with scanning nonlinear dielectric mi- 
croscopy, a width of only 0.5 nm in ultra-thin PZT films 
was determined fl2l ]. 

In the past ten years, piezoresponse force microscopy 
(PFM) has become a very common technique for domain 
imaging mainly due to its high lateral resolution with- 
out any need for specific sample preparation. In brief, 
for PFM a scanning force microscope is operated in con- 
tact mode with an alternating voltage applied to the 
tip. In ferroelectric samples this voltage causes thick- 
ness changes via the converse piezoelectric effect [l3| and 
therefore vibrations of the surface which lead to oscilla- 
tions of the cantilever that can be read out with a lock-in 
amplifier [II [lj]. 

In this contribution, we present a study of the influence 
of the tip radius on the PFM imaging of 180° ferroelec- 
tric domain walls in bulk single crystals. Therefore we 
determined the width W of domain walls when imaged 
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by PFM on LiNb03 crystals for tips of different radii. 
We also compared W for different types of crystals using 
one unique tip. Furthermore, a series of measurements 
was carried out with LiNb03 of different sample thick- 
nesses. Finally we present a simple analytical model that 
explains the observed dependencies of the domain wall W 
on the tip radius, the crystal thickness and type of sam- 
ple. 

The experiments were carried out with two different 
scanning force microscopes (Topomctrix Explorer from 
Veeco and SMENA from NT-MDT). Both systems were 
modified to allow application of voltages to the tip in or- 
der to enable PFM measurements. The alternating volt- 
age (f/tip = 10V PP ,30 — 60 kHz) was applied to the tip 
and the backside of the samples was grounded. To obtain 
reliable data for the width W we readout the A"-signal 
of the lock-in amplifier (in the following denoted as PFM 
signal). The recorded data were thus unaffected by the 
background inherent to PFM measurements [l6j which 
can lead to a presumed broadening of the width W [l7| • 

We used a series of different tips with radii varying 
from 10 — 90 nm, classified in Tab. HI All tips were made 
out of highly n-doped silicon and conductively coated 
with different materials. The spring constants ranging 
from 3 to 70 N/m. The tip with the smallest radius (r = 
10 nm) makes an exception since it was uncoatcd. Due 
to oxidization the outer few nanometers of the surface 
are modified to non-conductive Si02- As a consequence, 
the mechanical and the electrical tip do not coincide any 
more; the SiO^layer acts as a dielectric gap between tip 
and sample. 

The experiments for determining the dependence of 
the width W on the tip radius r were performed with 
periodically poled (period length A = 30 pm) congru- 
ently melting 500pm thick LiNbC>3 crystals (PPLN). For 
the measurements of W on the sample thickness t we 
used the same LiNbC>3 samples, mechanically thinned by 
polishing to the thickness wanted (15 — 1000 pm). The 
thinnest sample (0.9 pm) was a stoichiometric LiNb03 
crystal. Here the domains were generated with the help 
of the tip by applying a voltage of 20 V for lOmin. 
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TABLE I: Specifications of the different tips utiiized for the 
measurements. 



(a) 



Label 


Manufacturer 


Model 


Tip 


Coating 








radius [nm] 




A 


NT-MDT ° 


NSG11 


10 




B 


Veeco 6 


OSCM-PT 


15 


Ti-Pt 


C 


Veeco 6 


SCM-PIT 


20 


Pt-Ir 


D 


MikroMash c 


NSC35 


35 


Ti-Pt 


E 


MikroMash c 


NSC35 


50 


Cr-Au 


F 


NT-MDT a 


DCP11 


50-70 


diamond 


G 


MikroMash c 


NSC35 


90 


Co-Cr 



"www. ntmdt.ru 
'www.veeco.com 
c www. spmtips . com 



We also measured a series of samples different from 
LiNb03 as listed in Tab. [ill Those samples had thick- 
nesses of 0.5 — 2 mm and were either periodically poled 
(KTiOPC>4 and LiNbOa) or had arbitrary domain pat- 
terns (BaTi0 3 ,KNb0 3 ,LiTa0 3 and Sr .6iBa .39Nb 2 O 6 ). 

Calculating the spatial resolution achievable with PFM 
requires the exact electric field distribution underneath 
the tip and the electromechanical answer of the material. 
The latter is given by the dielectric constants as well as 
the elastic and piezoelectric tensors, respectively. This 
complex problem is most suitable for the finite element 
method (FEM), where all material constants can be in- 
cluded thus yielding quantitative results [HI ■ A detailed 
study on domain wall width imaged by PFM including 
FEM-calculations has recently been undergone [lj|. In 
this contribution, we propose a much simplified approach 
to the problem of lateral resolution in PFM. This model 
is not capable of giving the amplitude of the measured 
PFM signals because we performed some normalization 
in order to facilitate the calculations. The model can, 
however, give a quantitative prediction of the width W 
as a function of the tip radius r and the sample thick- 
ness t. 

In a first step the problem was simplified by approxi- 
mating the spherical apex of the tip with radius r by a 
point charge at the distance r above the sample surface. 
We further assumed the sample to be isotropic with an 
effective dielectric constant e e ff = \f£r £z, £z being the 
dielectric constant in z-dircction and e r the radial one 
perpendicular to z. The electric field distribution under- 
neath the tip is then given by (20| : 



E z (x,y,z) = — 

Seff 



z + r 



x 2 + y 2 + (z + r) 
z-r -It 



3/2 



(1) 



y 2 + (z - 2t - r) 



3/2 



Next we define E m!iX as the maximum electric field under- 
neath the tip apex inside the sample (at x = y = z = 0). 
fi max will be used to normalize the electric field strength. 




(b) 



active 
volume 



E.>E„ 



1 [ v" 






2 


J. J 



-4r -2r 2r 4r 
Tip radius r 




Positions; 



FIG. 1: (a) Electrical field distribution E z underneath a 
tip of radius r calculated with Eq. [1] for a sample of in- 
finite thickness. The white lines indicate the shells where 
the electric field has decayed to the indicated %-value of 
the maximum electric field S max present underneath the tip 
apex, (b) Schematics of the analytical model with the tip 
in the vicinity of a domain wall. The active volume where 
E z > Ecut ~ 0.05 -Emax underneath the tip covers both do- 
mains. The contributions of part I and III cancel out each 
other. The resulting surface deformation is thus determined 
by part II only. 



Furthermore, we scale all lengths with the tip radius r. 
The problem has thus become dimensionless. Figure[TJa) 
shows the electric field distribution E z inside a sample of 
infinite thickness. It can be seen, that at a depth of twice 
the tip radius r, the electric field has decayed to almost 
10% of its initial value E max . 

Figure [ljb) illustrates the analytical model used for 
the simulation of the lateral resolution of PFM depend- 
ing on the tip radius and the sample thickness. First we 
define an active volume inside of which the electric field 
has decreased to a value E cut which has to be identified 
later by fitting the experimental data. Beyond the active 
volume we set E z = 0. The resulting piczomechanical de- 
formation is calculated by integrating the contributions 
of the sample within the active volume: 



Az{x) 



133 z 



Ecu 



E z (x,y,z) 



dE y dE z 



(2) 



As can be seen in Fig. [TJb) the contributions of re- 
gion I and III to the emerging piczomechanical defor- 
mation cancel out each other. For the simulation of the 
PFM signal when scanning across a domain wall, the 
active volume was subdivided into approx. 10 9 cubic ele- 
ments. In order to determine E cnt the experimental data 
was fitted with the calculated slopes of the PFM signal 
across the domain wall. The best fits were obtained for 
-Ecut = 0.92 £ max for all tip radii. Note that the main as- 
sumption made for the simulation consists in the stiffness 
of the crystal, i.e., within a length of some microns the 
crystal parts can not deform independently. The strength 
of clamping in bulk samples has been demonstrated by 
other experiments [2lj | . 

In order to deduce the domain wall width W from the 
experimental data we normalized the PFM signals to an 
amplitude of 1 and fitted the data with a modified hy- 
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FIG. 2: Measured PFM signal line scan (•) across a 180° 
domain wall in LiNbOa recorded with a tip of type B (r = 
15 nm). The domain wall width W is determined fitting the 
data using Eq. [3j For comparison, a line scan is shown that 
was calculated with the theoretical model for a tip with a 
radius of 15 nm . 



FIG. 3: Measured domain wall width W as a function of the 
nominal tip radius r (Tab. [Jl. The straight line was calcu- 
lated using the analytical model presented in this contribu- 
tion. Note that despite its smaller radius, the uncoated tip A 
shows an inferior lateral resolution. 



pcrbolic tangent 

X(x) = Atanhf— ) + 5arctan( — ) (3) 

where A, B and w are used as free parameters. Note that 
this function is only used to determine the width W but 
has no direct physical meaning. We then applied a 25% — 
75% criterion on the scan lines of the PFM signal which 
corresponds to the full width at half maximum of the 
PFM amplitude if no PFM background is present [l7| . 
Figure [2] shows an example for a r = 15 nm tip on a 
LiNbC>3 sample. The measurement (dotted line) is fitted 
according to Eq. [3] (slim line) . For comparison the slope 
of the PFM signal calculated with the analytical model 
is also depicted (dashed line). As can be seen, the error 
determining W using Eq. [3] is minimal that is why we do 
not show any error bars in the subsequent graphs. 

Processing the data in the above described manner, we 
extracted the data for W as a function of the tip radius 
as shown in Fig. [3j The straight line results from our 
analytical model. The excellent agreement between mea- 
surement and model is striking and strongly sustains the 
model to give reasonable estimates on W. Implicit to our 
model, an atomically sharp domain wall is thus also sus- 
tained by the PFM measurements. Note that with tip B 
a width of only 17nm was measured, the smallest value 
recorded with PFM in bulk materials so far. This has to 
be compared with a recent publication, where a lowest 
limit for W of 65 nm in LiNbC>3 was estimated (22|. This 
value as well as our currently highest resolution is by no 
means a fundamental limit of the material itself but by 
the available tip sizes and imaging parameters. It is fur- 
thermore evidently seen that the non-coated tip A shows 
a substantially reduced lateral resolution (W ~ 50 nm) 
than pretended by its tip radius of only 10 nm. This, 
however, is exactly what can be expected from an sur- 
face oxide layer: the conductive part of the tip being 
at a distance of some nanometers from the sample sur- 
face, separated by the dielectric SiC>2-layer generates an 
less localized electrical field inside the crystal, thus the 
reduced spatial resolution. 



In a further series of experiments we determined W as 
a function of the sample thickness t varying the latter by 
three orders of magnitude from 0.9 - 1000 um (Fig. 2J). 
We used tips of type D with a nominal radius r = 35 nm. 
The accuracy and durability of the tips ware controlled 
by measuring a standard 500 um thick PPLN sample be- 
fore and after each data acquisition with a sample of 
modified thickness. From Fig. 2] no change of the width 
W within the thicknesses range of the samples could be 
observed. This, however, is consistent with our theoreti- 
cal model where the electric field distribution E z is found 
to be independent on the sample thickness t for t > 15r. 
In the case of tip D with a nominal radius of r = 35 nm 
the electric field E z is thus the same for samples with 
thicknesses > 500 nm. 

Finally we compared W for different crystals using tips 
of type B and D (Tab. [TTJ) . We again checked the reliabil- 
ity of the recorded data by always recording comparative 
measurements with a standard PPLN sample. All sam- 
ples show the same width W for a specifice tip of radius 
r = 15 nm or r = 35 nm within an error of ±lnm, al- 
though their dielectric anisotropy 7 = yj e z /e r differ as 
listed in Tab. fl] This can be understood if we calculate 
the electrical field c. g. with the method of image charges 
where we place the charge q at (0,0,— r). In contrast 
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■ 
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FIG. 4: Measured domain wall width W as a function of 
sample thickness t. The straight line is the result from the 
analytical model. 
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TABLE II: Domain wall width W measured for different sam- 
ples with tips of type B (r = 15 nm) and D (r = 35 nm). c- 
LiNbOs: congruently melting and s-LiNbOs: stoichiometric 
lithium niobate respectively, SBN: Sro.6iBao.39Nb2C>6. 



Sample 


Domain wall width 


Dielectric anisotropy 




W [nm] 


7 = y/sz/e r 




r — 15 nm / r = 35 nm 




BaTiOa 


19 / 46 


0.18 


KNbOa 


18 / 45 


0.34 


KTiOP0 4 


17 / 46 


0.86 


LiTa0 3 


18 / 45 


0.89 


c-LiNb0 3 


17 / 46 


0.58 


s-LiNb0 3 


17 / 45 


0.58 


Mg:LiNb0 3 


18 / 47 


0.58 


SBN 


18 / 48 


1.30 



to the case of an infinite isotropic half-plane, character- 
ized by a single dielectric constant, we need two image 
charges q' at (0, 0, r) and q at (0, 0, f) to account for the 
dielectric anisotropy given by e z and e r . The potential 
distribution above the surface is identical to the isotropic 
case but inside the crystal it is given by 



<p(z,y,z) = 



<1 



1 



(4) 



The only way to satisfy the boundary conditions at the 
surface (x, y, 0) is to set r = which leads to 



E z {x,y,z) 



2q 7 



X A _|_ yZ _|_ _|_ j, j 



3/2 



(5) 



From this equation we can clearly sec that the dielectric 
anisotropy does not affect the shape of the electrical field 
distribution but only the field strength. 

To summarize, we have analyzed the impact of the 
tip radius and the sample thickness on the domain wall 
width W observed with piezoresponse force microscopy. 
We introduced an analytical model which explains the 
experimental data: a linear dependency of W on the tip 
radius as well as no dependency of W on the sample 
thickness as long as the sample is thicker than 15-times 
the tip radius. The model assuming an infinite sharp do- 
main wall is perfectly consistent with the experimental 
data. Furthermore we carried out a series of measure- 
ments for different kind of single crystals. Even though 
their material parameters differ significantly all measured 
domain wall width were found to be independent of the 
specific sample. This is explained by the identical shape 
of the electrical field inside the sample which is the basis 
of our analytical model. 
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